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MODELING OF T.'UIN-FXUI GaAa GROWTH 


By 

John U. Hainbockel* 

* INTRODUCTION 

Hie Solid-on-SoUd «.SOS) raodel of cryacal gtowtiv (ref, ).) is repre- 
sented by a rectangular nrray of integers where each integer represents the 
number of adatoms in a coluun perpendicular to some reference frame. Tlie 
adatoms can represent atoms or molecules that are being stacked. Figure I 
illustrates the surface adatoms that are at the tops of their coluims. It 
is assumed that an adatom event of adsorption or desorption can only occur 
at the top of a column, 

We are concerned with constructing a model of crystal growth that takes 
into account the processes of nucleation on the growing surface as, well aa 
considering the processes of surface migration and dosorption of 
adatoms . 

In the SOS model the columns are constructed upon an M x M-square array 
by randomly placing adatoms upon the array and allowing these randomly 
deposited adatoms to either condense, evaporate, or migrate. Hie SOS model 
can be described as an array of interacting columns of varying integer 
heights. Tlie surface adatoms, being at the tops of coliunns, are allowed to 
migrate, remain stationary, or evaporate as is dictated by a set of rules 
which will be described presently.' 

Tlie term ’’epitaxy’' means ”an arrangement" and is used to denote the 
growth of one substance upon the crystal surface of a foreign substance. 

The term ’’autoepitaxy’’ is the oriented growth of a substance onto itself and 
"hetroepitaxy” is used for the growth of one material upon the surface of a 
different material. Obviously, hetroepitaxy becomes autoepitaxy after one 
layer of adatoms has been deposited over the growing surface. We uae the 
SOS method to simulate epitaxial growth of crystals, 

★Professor, Department of Mathematical Sciences, Old Dominion University, 
Norfolk, Virginia 23508, 
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The rules by which the columns of the SOS model interacted were gov- 
erned by the following ideas relating to the potential energy and potential 
energy changes associated with the adsorptioni migration, or desorption of 
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adatoros from an arl'/itrary row i and column j of an M x M array* Enct- 

giea aaaociated with an arbitrary site (i,j) wore defined aa follows; Uq * 
i\) (tij)“~tUe potential energy at a aite because of aurface bonding and crys- 
tal structure; (Ijo — the potential energy change at site (i,j) because of the 
deposition of an adatom (assumed the same for all sites); - Wj^(i ■ 1, 

8)— the potential energy changes at neighboring sites wlien an adatom is 
deposited at site (i^j); E(i|j)‘“‘the random surface energy associated with 
site (i,j) and time interval At; U(i,j) * Uo(i,j) + E(i,j)— the total ener- 
gy associated with aite (i,j) during the time interval At; Ug— the evapo- 
ration potential; and U„,--the migration potential. All of the above 
energies were measured in electron volts. 

We developed a Monte Carlo computer simulation of crystal growth (refs. 
2, 3, 4, and 5) by developing rules that determined the SOS kinetics of 
condensation evaporation or surface migration of adatoms. 'Iliese rules led 
to a cons.istent and physically reasonable description of the fundamentals 
associated with crystal growth. We first considered the adsorption of a 
thermally accommodated adatom onto the surface at some general site where 
the potential at this site was changed and, simultaneously, potential energy 
changes at all of the neighboring sites occurred. In Table I the potential 
energy changes are depicted by the mnemonic mask. Tlie center of this mask 
is placed over the site (i,j) to illustrate the changes to be made in the 
potential at the central site as well as the potential changes in the sur- 
rounding neighboring sites. 

The potential changes in the case of desorption of an adatom from the 
central aite are again depicted with the mask of Table 1, with the opposite 
signs on the potential changes. Tl\e case of surface migration was treated 
as a desorption from a site (i,j) followed by an adsorption at a nearest 
neighbor location, together with the correct potential mask changes associ- 
ated with each process. Ihe nearest neighbor migration site was determined 
by a random walk to one of the unoccupied nearest neighbor sites. 


Table 1. Potencial eueegy ehanges aaaociabed with central aice (ijj) 
ami neighbor altea due to depoaition of an adatom. 
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*nie Monte Carlo airaulation of crystal growth involved a random deposi- 
tion of thermally accoivunodated surface adatoms during a time interval At, 
11^080 deposited adatoiws changed the potential energies at the random surface 
sites under consideration. Tiie values assigned to the central potential 
change (J»q and neighboring potential changes -wi, i “ I, .i.» 8 
dictated the new potential energy values when an adatom was deposited or 
removed from a site, In this way each surface site had an energy barrier to 
translation or evaporation, represented by a potential well. We assuned 
that the thermally acconunodated adatoms had a surface energy distribution 
described by the Boltzmann distribution 

f(E) JL exp (-^), E > 0 
icr KT 


which has a mean energy of KT, 


During each time interval 
assign a random energy E(i,j) 

U(i,j) " U^Ciij) + E(i,j) 


At, the Boltjauann distribution 
to each of the surf ice adatoms. 


was used to 
We let 


( 2 ) 


denote the total energy possessed by a surface adatom at a site (i,j) 
during this time interval. This total energy is the sum of the potential 
energy due to the lattice structure and a random energy E from the 

Boltzmann distribution wi\ich characterizes the random surface energy. WIten 
U was leas than some material-dependent migration level the adatom 

remained stationary at the surface site. If U,„ < U < surface 
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migration by random walk was allowed to occur. If U was greater than the 
evaporation potential Ug, the adatom was removed from the site. 

Ihe rate of impingement of adatoms upon the surface was independent of 
the surface configuration. The rates associated with the evaporation and 
migration of adatoms depended upon the potential barriers Ug and Uni 
and also upon the values assigned to the potential changes and ~w£, 

(1 1, ..., 8). These later potential changes had to take into account the 

type of crystal structure and orientation of the growth we were trying to 
simulate with the SOS model. In Figure 2(a), for growth on the (100) face, 
we set up a correspondence between the central site, ti;q nearest neighbor 
potentials tj> , second nearest neighbor potentials ^nd the adatom 

potential changes for the mask in Table 1 (e.g,, wj (|>2 » W£ *■ (})i). In a 
similar manner we were able to set up the correspondences illustrated in 
Figure 2(b) and (c) for the (111) and (llO) orientations. In Table 2, we 

selected the relation between the neighbor potential? (fiQ, <^*1) <}>2» 'J'3 in 
such a way that when the first level of adatoms covered the surface, the 
potential distribution returned to its original value. By simply adding 
adatoms to the surface it was readily verified that the potential changes, 
assigned to the mask, had to adhere to the rules given in Table 2. In these 
rules, a negative sign denotes an attractive potential and after one 
complete layer of adatoms is deposited, the potential energy values at each 
site will return to their initial values. 


We let denote the change in the nearest neighbor potentials due to 

the addition of an adatom to the surface and let <j>2 > <j>3 denote the second 
and third nearest neighbor potential changes. We assumed that ((>2 * a2't'i 
and (|!3 "03(1)1 where 02, 03 are scale factors which are less than one. 

This allowed us to define the crystal orientation factor 5 as 
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which takes into account the different crystal orientations. We also defined 
the kink site potentials before *nd after and the capture of an 


.( 1 ) 


( 2 ) 


adatom as “ ^*^1* ^ks ** ^o 1 • (Note that " 2^*1'^.) 
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Table 2, Potential changes fov addition of an adatom to an arbitrary site. 
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THE SIMULATION MODEL AND PARAMETERS 


A flow chart of the simulation model is illustrated in Figure 3, The 
model is simple and presents an alternative viewpoint for the interaction of 
surface molecules: An assumed impingement rate dictates the number of 

adatoms arriving on the surface during a time interval At» Each of these 

t 

adatoms are added to the suri*ace at random sites and the potentials at each 
of these sites and neighboring sites are adjusted. If the At time 
interval is so small that no adatoms arrive on the surface, then every 
surface adatom can still be assigned a random energy from the Boltzmann 
distribution and surface interactions can be taken into account. We 
continued scanning the surface each At time interval until enough time 
accumulated for the addition of another adatom, 

The model allows for various .assumptions to be made about the 
interaction of potentials and assignment of potential values. We let 
» 0 denote the evaporation level, then. AUg « Ug - Uq represented the 
desorption energy AHevap, The activation energy for migration of adatoms 
in a flat surface was * U^ - Uq. The various potentials are 

illustrated in the Figure 4. The values assigned to U^ and U© 
greatly affected the model behavior. For example, the BoltiWann 
distribution is illustrated in Figure 5, where nominal values of AUm and 
AUg are illustrated. The number of surface adatoms with a statistical 
surface energy less than AUjjj is proportional to the area under the 
probability density curve which is given by N^ “ 1 - exp(-AU^/KT) , 

The number of adatoms that escaped from the surface is proportional to 
N « exp(-AU /KT) and the number of adatoms that migrated is proportional 

to the area N » I “ N. - N . Letting a =AU /AU , Figure 6 was con- 
m £ e ® m e' ® 

structed which illustrates the migration effect as a decreases. 

The values of (})q , ({ij , 't>3 which denote the potential energy 

changes at a central site (i,j) and nearest neighbor sites can be different 
for the substrate and the growing material. For the substrate material we 
could use the depth of the surface potentials and migration levels to 
stimulate a variety of surface morphologies. In this model we envisioned a 
flat substrate as a periodic lattice structure 20-X-20 square where each 
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lattice Is a potential well* The substrate can vary from flat to rough and 
the potentials adjusted to reflect various surface preparations. For an 
ideally flat substrate we assumed that the depths of the potential wells 
were uniform, given by Uos • After one layer of growing material 
covered the surface, the potentials at each site were assumed to convert to 
the autoepitaxy potentials IVy* In order to moke this transition we 


assumed that 


8 

- I w. + (U 

O I 0 


Uo,)r.. where 



is zero if the 


height h|j at position (ij) is greater than or equal to one and 
Tj[j is one in the case where h-j[j is zero. Thus, if an adatom was 
deposited at a first layer site (i,j) we adjusted the potential at this site 
by the relation by *• in addition to the mask potential changes as 

this produced the desired change that hetroepitaxy produces in the potential 
at the surface site. 


Nucleation on the substrate was controlled by the values assigned to 

AU AU , and <|> , For large values of AU there were deep potential 
e, m o o jjj 

wells that captured all thermally accommodated adatoms. For small AU,j, 

values there was an increase in surface migration and a decrease in the 

number of adatoms that remained localized, This increased the probability 

of an adatom combining with other adatoms to form a critical cluster. Then 

growth was characterized by the lateral motion of adatoms and their addition 

to the steps of clusters that produced the lateral growth. 


Various potentials were proposed for the addition of an adatom to the 
surface (refs. 6 and 7): 

Buckingham Potential E ■ He”®^ 1 


Modified Buckingham Potential E " 

Lennard- Jones (Mie Potential) E « 

Morse Potential E ■ ” ^o^ 

BornrMayer Potential E ^ Ae”®*^ 







rm 
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These potentials reflect the vertical effect of potential change. For 
the lateral Interaction between potentials and resultant change?* (ref. b), 
we find: 

Kiselev Potential B " E - NCf(r) 

0 

where is interaction at zero coverage, G is dispersion constant, N 
is the number of nearest neighbors at half a monolayer coverage and r is 
mean distance between molecules. 

Output from the computer program can be graphic as illustrated in the 
Figure 7 or quantitative. 


QUANTITATIVE MEASURES OF CRVSTAI GROWTH AM) PARAMETERS OF MODEL 


Measures o£ Cryglu^i Growth 


P a raroe t era of Model 


1, Growth rate of crystal 

2, Critical clusters 

a. sizQ 

b. shape 

c. density 

Vs, time or deposition rate 

3, Surface diffusion (mobility) 

4, Condensation rate 
5# Evapoi'ation rate 

6. Rate of nucleation 

7. Other characteristics 


U 

2 , 

3, 


4. 

5, 
6# 
1 , 
8 . 

9. 


10 , 


11 , 


Deposition rate of adatoms 
Potential changes <{)qi ^3 

Mean Uq and standard deviation 
Oq associated with normal distri- 
bution N(Uo» <Jo) of surface 
potentials (initially Uq " 0) 

Traps in Surface 

Temperature of substrate 

Number of migration scans (time At) 

Crystal orientation 

Substrate and growing potentials can 

be different 


Mean Urn and standard deviation 
Ora ot migration levels associated 
with normal distribution NiU^, 


Ora) . 

Initial substrate geometry and 
potentials 

Assumptions in regard to retention 
of incident energy 


IB Is * 

*^random “Boltsraann ‘^Retention 
(Surface (Kinetic 

Energy) Energy of 

Incident 
Ad at om) 


NOMINAL VALUES FOR POTENTIAL ENERGIES FOR 
GERMANIUM IN EV (1 EV - 23 KGAL/MOLE)** 
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ev 

Ge 
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Ge 
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Ge 
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See ref. 9, 
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Figure 1. SOS model for crystal growth. 
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Figure 2. FCC model and potential chang( 
crystal orientations. 
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Figure 3. Plow chart of simulation model. 
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Figure 7. Graphic display of crystal growth (100) 
orientation. 
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